Calibration and Characterization of the IceCube 
Photomultiplier Tube 



R. Abbasi 2 , Y. Abdou\ T. Abu-Zayyad af , J. Adams , J. A. Aguilar 2 , 
M. Ahlers ae , K. Andeen 2 , J. Auffenberg al , X. Bai ad , M. Baker 2 , 
S. W. Barwick v , R. Bay g , J. L. Bazo Alba am , K. Beattie h , J. J. Beatty°- r , 
S. Bechet 1 , J. K. Becker-*, K.-H. Becker al , M. L. Benabderrahmane am , 
J. Berdermann am , P. Berghaus 2 , D. Berley p , E. Bernardini am , D. Bertrand 1 , 

D. Z. Besson x , M. Bissok a , E. Blaufuss p , D. J. Boersma a , C. Bohm ag , 
O. Botner aj , L. Bradley ai , J. Braun 2 , S. Buitink h , M. Carson*, D. Chirkin 2 , 
B. Christy p , J. Clem ad , S. Cohen w , C. Colnard", D. F. Cowen ai > ah , 
M. V. D'Agostino g , M. Danninger ag , C. De Clercq m , L. Demirors w , 
O. Depaepe 111 , F. Descamps*, P. Desiati 2 , G. de Vries-Uiterweerd* , 
T. DeYoung ai , J. C. Diaz-Velez 2 , J. Dreyer s 'j, J. P. Dumm 2 , 
M. R. Duvoort ak , R. Ehrlich p , J. Eisch 2 , R. W. Ellsworth 15 , O. Engdegard aj , 

S. Euler a , P. A. Evenson ad , O. Fadiran d , A. R. Fazely f , T. Feusels*, 
K. Filimonov g , C. Finley ag , M. M. Foerster ai , B. D. Fox ai , A. Franckowiak 1 , 
R. Franke am , T. K. Gaisser ad , J. Gallagher 7 , R. Ganugapati 2 , M. Geisler a , 
L. Gerhardt h ' g , L. Gladstone 2 , A. Goldschmidt h , J. A. Goodman p , 
D. Grant ai , T. Griesel ab , A. Grofi°' u , S. Grullon 2 , R. M. Gunasingha f , 
M. Gurtner al , C. Ha ai , A. Hallgren a j, F. Halzen 2 , K. Han°, K. Hanson 2 , 
Y. Hasegawa 11 , J. Haugen 2 , K. Helbing al , P. Herquet ac , S. Hickford , 
G. C. Hill 2 , K. D. Hoffman p , A. Homeier 1 , K. Hoshina 2 , D. Hubert 111 , 
W. Huelsnitz p , J.-P. HiilB a , P. O. Hulth ag , K. Hultqvist ag , S. Hussain ad , 
R. L. Imlay f , M. Inaba 11 , A. Ishihara 11 , J. Jacobsen 2 , G. S. Japaridze d , 
H. Johansson ag , J. M. Joseph 11 , K.-H. Kampert al , A. Kappes 2 ' 1 , T. Karg al , 



* Corresponding author 
Email addresses: mase@hepburn. s . chiba-u. ac . jp (K. Mase), 
chwendt@icecube.wisc.edu (C. Wendt), syoshida@hepburn. s . chiba-u. ac . jp 
(S. Yoshida) 

1 affiliated with Universitat Erlangen-Nurnberg, Physikalisches Institut, D-91058, Er- 
langen, Germany 

2 on leave of absence from Universita di Bari and Sezione INFN, Dipartimento di Fisica, 
1-70126, Bari, Italy 

3 afhliated with School of Pure and Applied Natural Sciences, Kalmar University, S- 
39182 Kalmar, Sweden 

4 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 



Preprint submitted to Nuclear Instruments and Methods A 



February 11, 2010 



A. Karle z , J. L. Kelley z , N. Kemming 1 , P. Kenny x , J. Kiryluk h ' g , 
F. Kislat am , N. Kitamura 2 , S. R. Klein h ' g , S. Knops a , G. Kohnen ac , 
H. Kolanoski', L. K6pke ab , D. J. Koskinen ai , M. Kowalski k , T. Kowarik ab , 
M. Krasberg 2 , T. Krings a , G. Kroll ab , K. Kuehn q , T. Kuwabara ad , 
M. Labare 1 , S. Lafebre ai , K. Laihem a , H. Landsman 2 , R. Lauer am , 
A. Laundrie 2 , R. Lehmann 1 , D. Lennarz a , J. Liinemann ab , J. Madsen af , 
P. Majumdar am , R. Maruyama 2 , K. Mase 11 '*, H. S. Matis h , M. Matusik al , 
K. Meagher p , M. Merck 2 , P. Meszaros ah ' ai , T. Meures a , E. Middell am , 
N. Milke 8 , H. Miyamoto 11 , T. Montaruli 2 ' 2 , R. Morse 2 , S. M. Movit ah , 
R. Nahnhauer am , J. W. Nam v , U. Naumann al , P. NieBen ad , D. R. Nygren h , 
S. Odrowski u , A. 01ivas p , M. 01ivo aj ' j , M. Ono 11 , S. Panknim, L. Paul a , 

C. Perez de los Heros aj , J. Petrovic 1 , A. Piegsa ab , D. Pieloth 8 , A. C. Pohl aj ' 3 , 
R. Porrata g , J. Posselt al , P. B. Price g , M. Prikockis ai , G. T. Przybylski h , 
K. Rawlins c , P. Redl p , E. Resconi u , W. Rhode 8 , M. Ribordy w , A. Rizzo m , 

P. Robl aa , J. P. Rodrigues 2 , P. Roth?, F. Rothmaier ab , C. Rott q , 
C. Roucelle u , D. Rutledge ai , B. Ruzybayev ad , D. Ryckbosch*, 
H.-G. Sander ab , P. Sandstrom 2 , S. Sarkar ac , K. Schatto ab , S. Schlenstedt am , 
T. Schmidt p , D. Schneider 2 , A. Schukraft a , A. Schultes al , O. Schulz u , 
M. Schunck a , D. Seckel ad , B. Semburg al , S. H. Seo ag , Y. Sestayo u , 
S. Seunarine , A. Silvestri v , A. Slipak ai , G. M. Spiczak af , C. Spiering am , 
M. Stamatikos q ' 4 , T. Stanev ad , G. Stephens ai , T. Stezelberger h , 
R. G. Stokstad 11 , S. Stoyanov ad , E. A. Strahler 111 , T. Straszheim p , 
G. W. Sullivan p , Q. Swillens 1 , I. Taboada c , A. Tamburro af , O. Tarasova am , 
A. Tepe c , S. Ter-Antonyan f , C. Terranova w , S. Tilav ad , P. A. Toale ai , 

D. Tosi am , D. Turcan p , N. van Eijndhoven 111 , J. Vandenbroucke g , 
A. Van Overloop*, J. van Santen 1 , B. Voigt am , D. Wahl aa , C. Walck ag , 
T. Waldenmaier 1 , M. WallrafP, M. Walter am , C. Wendt 2 -*, S. WesterhofP, 
N. Whitehorn 2 , K. Wiebe ab , C. H. Wiebusch a , G. Wikstr6m ag , 

D. R. Williams b , R. Wischnewski am , H. Wissing p , K. Woschnagg g , C. Xu ad , 

X. W. Xu f , G. Yodh v , S. Yoshida 11 -*, P. Zarzhitsky b 

a III. Physikalisches Institut, RWTH Aachen University, D-52056 Aachen, Germany 
b Dept. of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA 
c Dept. of Physics and Astronomy, University of Alaska Anchorage, 3211 Providence Dr., 

Anchorage, AK 99508, USA 
d CTSPS, Clark- Atlanta University, Atlanta, GA 30314, USA 
e School of Physics and Center for Relativistic Astrophysics, Georgia Institute of 
Technology, Atlanta, GA 30332. USA 
Wept, of Physics, Southern University, Baton Rouge, LA 10813, USA 



2 



9 Dept. of Physics, University of California, Berkeley, CA 94720, USA 
h Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 
l Institut fiir Physik, Humboldt- Universitat zu Berlin, D- 12489 Berlin, Germany 
i Fakultat fur Physik & Astronomie, Ruhr- Universitat Bochum, D-44 780 Bochum, 

Germany 

k Physikalisches Institut, Universitat Bonn, Nussallee 12, D-53115 Bonn, Germany 
l Universite Libre de Bruxelles, Science Faculty CP230, B-1050 Brussels, Belgium 
m Vrije Universiteit Brussel, Dienst ELEM, B-1050 Brussels, Belgium 
n Dept. of Physics, Chiba University, Chiba 263-8522, Japan 
"Dept. of Physics and Astronomy, University of Canterbury, Private Bag 4800, 

Christchurch, New Zealand 
p Dept. of Physics, University of Maryland, College Park, MD 20742, USA 
q Dept. of Physics and Center for Cosmology and Astro- Particle Physics, Ohio State 
University, Columbus, OH 43210, USA 
r Dept. of Astronomy, Ohio State University, Columbus, OH 43210, USA 
s Dept. of Physics, TU Dortmund University, D-44221 Dortmund, Germany 
t Dept. of Subatomic and Radiation Physics, University of Gent, B-9000 Gent, Belgium 
u Max- Planck- Institut fur Kernphysik, D-69177 Heidelberg, Germany 
""Dept. of Physics and Astronomy, University of California, Irvine, CA 92697, USA 
w Laboratory for High Energy Physics, Ecole Polytechnique Federate, CH-1015 Lausanne, 

Switzerland 

x Dept. of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA 
v Dept. of Astronomy, University of Wisconsin, Madison, WI 53706, USA 
z Dept. of Physics, University of Wisconsin, Madison, WI 53706, USA 
aa Physical Sciences Laboratory, University of Wisconsin, Madison, WI 53706, USA 
ab Institute of Physics, University of Mainz, Staudinger Weg 7, D-55099 Mainz, Germany 
ac Universite de Mons, 7000 Mons, Belgium 
ad Bartol Research Institute and Department of Physics and Astronomy, University of 

Delaware, Newark, DE 19716, USA 
ae Dept. of Physics, University of Oxford, 1 Keble Road, Oxford 0X1 3NP, UK 
af Dept. of Physics, University of Wisconsin, River Falls, WI 54022, USA 
ag Oskar Klein Centre and Dept. of Physics, Stockholm University, SE-10691 Stockholm, 

Sweden 

ah Dept. of Astronomy and Astrophysics, Pennsylvania State University, University Park, 

PA 16802, USA 

al Dept. of Physics, Pennsylvania State University, University Park, PA 16802, USA 
a Wept. of Physics and Astronomy, Uppsala University, Box 516, S-75120 Uppsala, 

Sweden 

ak Dept. of Physics and Astronomy, Utrecht University /S RON, NL-3584 CC Utrecht, The 

Netherlands 

al Dept. of Physics, University of Wuppertal, D-42119 Wuppertal, Germany 
am DESY, D- 15735 Zeuthen, Germany 



3 



Abstract 



Over 5,000 PMTs are being deployed at the South Pole to compose the 
IceCube neutrino observatory. Many are placed deep in the ice to detect 
Cherenkov light emitted by the products of high-energy neutrino interactions, 
and others are frozen into tanks on the surface to detect particles from atmo- 
spheric cosmic ray showers. IceCube is using the 10-inch diameter R7081-02 
made by Hamamatsu Photonics. This paper describes the laboratory char- 
acterization and calibration of these PMTs before deployment. PMTs were 
illuminated with pulses ranging from single photons to saturation level. Pa- 
rameterizations are given for the single photoelectron charge spectrum and 
the saturation behavior. Time resolution, late pulses and afterpulses are 
characterized. Because the PMTs are relatively large, the cathode sensitiv- 
ity uniformity was measured. The absolute photon detection efficiency was 
calibrated using Rayleigh-scattered photons from a nitrogen laser. Measured 
characteristics are discussed in the context of their relevance to IceCube event 
reconstruction and simulation efforts. 

Keywords: PMT, neutrino, cosmic rays, ice, Cherenkov 
PACS: 85.60.Ha, 96.40.Z 



1. Introduction 

IceCube [H, 0] is a kilometer-scale high energy neutrino telescope currently 
under construction at the geographic South Pole. A primary goal is to detect 
high energy neutrinos from astrophysical sources, helping to elucidate the 
mechanisms for production of high energy cosmic rays 



IceCube uses the 2800 m thick glacial ice sheet as a Cherenkov radia- 
tor for charged particles, for example those created when cosmic neutrinos 
collide with subatomic particles in the ice or nearby rock. Neutrino interac- 
tions can create high energy muons, electrons or tau particles, which must 
be distinguished from downgoing background muons based on the pattern 
of light emitted. The Cherenkov light from these particles is detected by 
an embedded array of Digital Optical Modules (DOMs), each of which in- 
corporates a 10" diameter R7081-02 photomultiplier tube (PMT) made by 
Hamamatsu Photonics. The DOMs transmit time-stamped digitized PMT 
signal waveforms to computers at the surface. 

The finished array will consist of 4800 DOMs at depths of 1450-2450 m, 
deployed at 17 m intervals along 80 vertical cables, which in turn are arranged 
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in a triangular lattice with a horizontal spacing of approximately 125 m. An 
additional 320 DOMs will be frozen into 1.8 m diameter ice tanks located at 
the surface to form the IceTop array, which is designed for detection of cosmic 
ray air showers. The geometrical cross sectional area will be ~1 km 2 and the 
volume of ice encompassed will be ~lkm 3 . Another 360 DOMs will be 
deployed in a more compact geometry ("Deep Core" fi|) using PMTs almost 
identical to those described here but with a higher efficiency photocathode. 

In this paper we describe measurements characterizing and calibrating 
IceCube PMTs, and discuss their relevance to detector performance and 
event reconstruction. First we describe the signals of interest in Section [2j 
Section [3] briefly describes the DOMs in which IceCube PMTs are deployed. 
Section @] describes selection and basic features of the PMT, including the 
dark noise rate. Section [5] presents the design of the HV divider circuit. Sec- 
tions [6]-[9] discuss characteristics of the PMT in the photon counting regime, 
starting with single photon waveforms and charge distributions. Time resolu- 
tion is studied with a pulsed laser system. Uniformity of the photon detection 
response on the photocathode area is measured by scanning the entire cath- 
ode surface with a UV LED. Absolute efficiency calibration of the IceCube 
PMTs is carried out using Rayleigh-scattered light from a calibrated laser 
beam. Sections [T0T[TT1 describe response to bright pulses of light, including 
saturation behavior and afterpulse characteristics. 



2. Characteristics of optical signals in IceCube 

We begin by summarizing what the PMTs are supposed to detect, namely 
the optical signals generated by neutrinos in IceCube [H, 0]. Of particular 
relevance are the amplitudes and widths of the pulses, requirements on time 
resolution, and how the pulses are used to reconstruct physics events or reject 
backgrounds. 

In detection of a high energy by IceCube, the neutrino interaction 
creates a muon that traverses kilometers of ice and generates Cherenkov 
light along its path. Above 1 TeV, the muon loses energy stochastically to 
produce multiple showers of secondary particles, resulting in an overall light 
yield proportional to the muon energy i|, @] . Most light is emitted near the 
Cherenkov angle, which is 41° away from the track direction. The arrival 
times of detected photons depend on the position of each DOM relative 
to the muon's path. For close DOMs, most photons arrive in a pulse less 
than 50 ns wide, in which the earliest photons have traveled straight from 
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the muon track without scattering. Significant scattering accumulates along 
photon trajectories with a characteristic length scale of about 25 m Q, so 
light pulses lengthen with distance and reach 1 /j,s (FWHM) for DOMs 160 m 
away from a muon track. Depending on primary energy and distance from 
the track, each PMT can see single photons or pulses ranging up to thousands 
of photons. 

Event reconstruction [l], 0,0 builds on principles established for the pre- 
decessor array, AMANDA @, LLOfl . The observed PMT waveforms from indi- 
vidual DOMs are correlated and built into events, which are fitted to physics 
hypotheses using the maximum likelihood method. Each fit has access to 
the complete pattern of light amplitude and timing seen by the DOMs, and 
accounts for the DOMs' time response and optical sensitivity as well as time 
dispersion and optical attenuation introduced by the ice. The fit gives the 
direction and energy of the muon, which in turn characterizes its parent 
neutrino. 

The observed light pattern is also used to distinguish the rare neutrino 
events from the large background of muons created in cosmic ray air showers, 
which are 10 6 times more numerous. For this the reconstructed direction is 
key, because neutrinos can come from any direction, even up from below, 
but the background muons are downgoing. A small fraction of background 
events can be misreconstructed in direction, thereby appearing to come from 
neutrinos, but the pattern of detected light will generally be a poor match 
compared to expectations for a properly reconstructed track. Misreconstruc- 
tion can be aggravated by additional muons from the same shower or other 
coincident showers. This separation between signal and background is ac- 
complished by evaluating relative probabilities on an event-by-event basis, 
and is aided by good time and amplitude resolution as well as by low PMT 
noise rates. 

Similar principles apply to other types of high energy neutrino interac- 
tion. Instead of a muon, an electron can be created that loses its energy in 



a few-meter- long particle shower ; on the scale of IceCube, such a shower 
appears almost like a point source of Cherenkov light. For a sufficiently ener- 
getic neutrino, the light can be detected hundreds of meters away, and nearby 
DOMs can see enough light to drive their PMTs into the nonlinear saturation 
regime. Therefore proper modeling of saturation behavior is needed for good 
reconstruction and background rejection. 

Design studies j^j for important physics goals have shown that sufficient 
reconstruction quality is achieved for a PMT timing resolution of 5 ns, low- 



6 



temperature noise rate below 500 Hz, and effective dynamic range of 200 
photoelectrons per 15 ns. 

In the case of lower energy (MeV) neutrinos from supernovae, IceCube 
cannot resolve individual interactions. Instead, supernovae would be detected 
as a momentary increase in the collective photon counting rate for the whole 
array, corresponding to a large number of neutrino interactions within a few 
seconds. The dark noise rate of the PMTs is particularly important here, 
because it dictates the statistical significance of any excess count rate. 

IceTop uses DOMs identical to those in the deep ice. Here the signals 



arise from muons, electrons and gamma rays in cosmic ray air showers [12 



These particles deposit energy in the ice tanks housing the DOMs, resulting 
in light pulses up to several hundred nanoseconds long. The arrival times 
and amplitudes in the surface array are then used to reconstruct the shower 
core position, direction, and energy. An overall timing resolution of 10 ns 
provides pointing accuracy of about 1 degree. The PMT pulses range from 
single photoelectrons at the periphery of showers to 10 5 photoelectrons for a 
1 EeV shower that strikes within the array. To achieve the implied dynamic 
range, each tank contains two DOMs operating at gains differing by a factor 
50. 



3. The IceCube optical detector: DOM 

The Digital Optical Module is the fundamental element for both opti- 
cal detection and data acquisition in IceCube 0, 0, 14]. It contains a 10" 



diameter PMT supported by coupling gel, the PMT high voltage generator 
and divider circuits, an LED flasher board used for calibration of the ar- 
ray geometry and study of ice properties, and the DOM mainboard which 



contains the analog and digital signal processing electronics [14J. The PMT 
is surrounded by a mu-metal grid to shield it from the terrestrial magnetic 
field and improve the PMT performance. All systems are housed within a 
pressure sphere made of 0.5" thick glass, capable of withstanding pressures 
to 70 MPa. The glass and gel set the short wavelength cutoff of the DOM at 
about 350 nm, where the PMT by itself is still relatively efficient. 

Strings of DOMs are deployed into water columns that have been melted 
by a hot-water drill. After refreezing, DOMs are optically well coupled to 
the surrounding glacial ice. Signal and power connections between the DOMs 
and the surface are provided by copper twisted-pair wires bundled together to 
form the main cables. PMT signals are digitized on the mainboard, buffered 
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Table 1: Hamamatsu specifications for the R7081-02 PMT (typical) [lj|. 



Spectral response 300 to 650 nm 

Quantum efficiency at 390 nm 25% 

Supply voltage for gain 10 7 1500 V 

Dark rate at -40 °C 500 Hz 

Transit time spread 3.2 ns 
Peak to valley ratio for single photons 2.5 

Pulse linearity at 2% deviation 70 mA 



in memory, and sent to the surface upon command of surface readout pro- 
cessors. 



4. PMT selection and dark noise rate 

A number of large-area PMTs are commercially available and have been 
used successfully to instrument large volumes in other experiments. IceCube 
selected the R7081-02 made by Hamamatsu Photonics, emphasizing the cri- 
teria of low dark noise and good time and charge resolution for single photons. 
Some manufacturer's specifications are shown in Table [TJ and more detailed 
measurements are described in the following. 

The nominal gain of 10 7 was chosen to give single photon pulses around 
8mV, which is well above the digitizer precision and other electronic noise 
levels (both ~0.1mV). Aging was not a concern for gain selection, since at the 
expected noise rates, the corresponding total charge delivered by each deep- 
ice PMT will be less than 1C after 20 years (or 100 C for Ice Top). Tubes 
with 10 and 12 stages were evaluated, with the 10 stage options showing a 
better peak-to-valley ratio at this gain. Lower gains of 5 x 10 6 and 10 5 were 
chosen for Ice Top PMTs because air shower pulses generally comprise many 
photon detections. 

The R7081-02 has 10 linear focused dynode stages and achieved the nom- 
inal gain of 10 7 at about 1300 V in our tests with the recommended divider 
ratios (observed range 1050-1600 V for 3744 PMTs). Its 10" diameter photo- 
cathode is composed of the standard bialkali material (Sb-Rb-Cs, Sb-K-Cs) 
with a peak quantum efficiency of approximately 25% at 390 nm. With a 



borosilicate glass envelope, the spectral response [16| is a good match to 
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the spectrum of Cherenkov light after propagation through ice [7|, especially 
considering the 350 nm cutoff of the housing. 
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Figure 1: Dark count rate versus temperature for a sample IceCube PMT covered with 
black tape (see text). Rates were recorded after a settling time of 12 hours of dark 
operation at gain 10 7 and discriminator threshold 0.25 times the single photoelectron 
peak. An artificial deadtime circuit rejected additional hits within 6 [is of each count, 
including about half of all afterpulses (Section ITT]) . The contribution from cosmic rays 
(<5 Hz) has not been subtracted. 



In response to IceCube requirements, the supplied R7081-02 units were 
manufactured with a custom low radioactivity glass. The resulting dark 
count rate at low temperatures is close to 300 Hz in the — 40 °C to — 20 °C 
range of greatest interest for IceCube (Fig. [TJ). The higher room tempera- 
ture rate can be attributed mostly to cathode thermionic emission, which is 
suppressed at low temperature. The low temperature rate is believed to be 
dominated by radioactive decays plus scintillation in the PMT glass enve- 
lope, and shows a rise with decreasing temperature similar to that reported 



in other studies [17]. The association of this rate with decays is supported 
by time correlations observed on scales up to 1 msec, as can result from 
delayed particle capture or de-excitation of states created by decays. It is 
also supported by the observed effect of taping: for these measurements, 
the entire outside surface of the PMT glass was covered in black vinyl tape, 



9 



pulled tightly against the glass to avoid bubbles. The taping is observed to 
reduce the low-temperature noise rate by about half. The reduction is at- 
tributed to absorption of outward going decay photons, which can otherwise 
be channeled to the photocathode via internal reflection. The taped result is 
appropriate for PMTs installed in IceCube DOMs because they are optically 
coupled to gel (then glass and ice) where the refractive index matches better 
than it does for air. 

The low dark rate allows IceCube to record all events that satisfy simple 
multiplicity conditions, and is particularly important for observation of any 
galactic supernova event. Such a supernova could yield about 10 6 excess 



photon counts in IceCube over a few seconds [18|. The single-PMT dark 
rate, multiplied by the number of PMTs, contributes a background rate of 
1.5 x 10 6 Hz, with a similar contribution from decays in the DOMs' glass 
housing. The excess from a supernova would be easily observed above this 
background, even including the details of its time structure. 

On the other hand, a high energy neutrino event creates optical pulses 
distributed over 3 /xs, with most information contained within a time window 
less than 300 ns wide in each DOM. Because this window is so short, the low 
PMT dark rate implies that only 1% of muons would be accompanied by a 
relevant noise count among the 100 DOMs closest to the track, and many 
of these DOMs detect multiple signal photons. Therefore the degradation 
of reconstruction and background rejection is very small. The dark noise 
rate has even less effect for IceTop DOMs, due to higher thresholds and 
coincidence requirements. 

5. High voltage divider circuit 

The relative dynode voltage ratios for R7081-02 have been optimized by 
Hamamatsu to achieve a maximum collection efficiency while achieving 10 7 
gain between 1050 V and 1600 V. Our high voltage subsystem design fulfills 
the additional requirements of low power consumption, long-term reliability, 
and sustained response to very bright pulses lasting up to a microsecond. 

The dynode voltages are provided by a passive resistive divider with a 
total resistance of 130 MQ (Fig. [5]). The rather high total resistance is chosen 
to minimize power consumption, which is an important economic consider- 
ation for operations at the South Pole. A custom, compact, high voltage 
generator 1^] that is both low power (< 300 mW) and low noise (< 1 mV 



ripple, peak-to-peak) is used in conjunction with the passive divider. 
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Figure 2: Schematic of the passive HV divider circuit, shown with coupling to the front-end 
amplifier on the digitizer board. 



Capacitors are placed across the last six dynode intervals and between 
the last dynode and anode. These capacitors help sustain the PMT output 
for large pulses of up to 10 6 photoelectrons (p.e.). Even for illumination 
in the PMT saturation region (~200 p.e./ns, see Section [TO]) , the transient 
gain loss after a 1/is pulse (2 x 10 5 p.e.) is observed to be less than 1%. 
A detailed simulation |20| indicates that such a pulse could arise from a 
50PeV electron shower 100 m away from the PMT, which would then be 
faithfully recorded. Pulses up to five times this large (10 6 p.e.) still result 
in less than 20% transient gain loss, so while the primary pulse would be 
completely saturated, afterpulse amplitude could be used to estimate the 
total illumination (Section [TTj) . Finally above 2 x 10 7 p.e. the gain loss rises 
rapidly above 50%. For all these transient gain losses, recovery occurs within 
the RC time constants of order 1 sec. 

Low-inductance resistors (100 Q, R13 through R15) are used to dampen 
ringing that arises from coupling of the larger dynode filter capacitors with 
parasitic inductance in the dynode leads and printed circuit traces. This 
ringing could otherwise be a nuisance when reconstructing a single PMT 
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output waveform as a series of photon hits. 

The IceCube PMTs are operated with their cathode at ground potential. 
Therefore the high voltage anode is AC coupled to the front-end amplifiers. 
For the AC coupling, we use a custom bifilar-wound 1:1 toroidal transformer 
rather than a DC blocking capacitor. High voltage reliability is achieved 
in the transformer winding using wire with insulation rated for over 5kV 
DC. The resulting stray capacitance from anode to front-end amplifier is 
only 30 pF, which limits the stored energy which might damage the analog 
front-end if sparking should occur in the HV system. In contrast, a cou- 
pling capacitor large enough to meet the signal droop specification would 
be at least 1000 times larger than the stray capacitance of the transformer. 
The transformer topology also reduces noise by avoiding noisy high voltage 
ceramic signal coupling capacitors and by breaking a ground loop path in- 
volving the HV power supply. The ferrite and resistors in series with the 
HV supply further reduce coupling of high frequency noise to the front end 
input. 

The transformer coupling delivers good signal fidelity for single-photoelectron 
(SPE) waveforms with risetimes of a few ns, while transmitting wide pulses 
exceeding 1 /is with less than 10% droop and undershoot. The custom design 
uses 18 bifilar turns on a ferrite (Magnetics Type H) toroid core, providing 
roughly flat coupling from 8 kHz to over 100 MHz at operating temperatures 
down to — 40 °C. The self-resonant frequency is above 150 MHz. The low 
operation temperature presented a challenge because the permeability of the 
transformer core decreases rapidly with temperature, leading to a shorter 
droop time constant. Although most DOMs have a time constant around 
15 /is at ambient temperatures near —30 °C, 1200 DOMs were built using an 
older transformer design yielding a time constant around 1.5 /xs at — 30 °C. 
The improved performance of the new design was achieved with a larger core 
and more turns, at the expense of a slightly wider SPE pulse shape (Fig. [3]). 
The two designs are deployed intermixed. 

The droop and undershoot are relevant for the /is long trains of photon 
pulses expected in DOMs over 100m away from high energy events, such 
as lOTeV electron showers or 500 TeV muon tracks. The small remaining 
effects are corrected by a software digital filter as a first step in event recon- 
struction, based on individual time constants for each DOM. The residual 
error is typically less than 1% of the pulse amplitude (except for pulses with 
peak or undershoot outside the ADC dynamic range, which is limited after 
400ns 0). 
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The divider circuit is constructed on a 10 cm-diameter printed circuit 
board which is directly solder-mounted to the PMT. All components (except 
R13-R16) are through-hole mount type, selected with a voltage derating 
factor of two or greater (typically four) to ensure long-term reliability. Strict 
voltage and voltage gradient rules are applied to the board layout. 

Coaxial cables are used for the connections to the high voltage generator 
and the front end amplifier on the digitizer board. The effective load for 
anode output pulses is 50 Q (43 Q for the older transformer design), which 
includes a back-termination resistor on the primary side of the transformer 
(R16), the transformer AC response, and the input impedance of the ampli- 
fier. 

6. Single photoelectron waveform and charge 

The SPE waveform shape and charge probability distribution are impor- 
tant for event reconstruction. The DOM waveform digitizers are triggered 
when the signal reaches about 0.25 times the typical SPE peak amplitude, 
after which the PMT output waveform is digitized for up to 6.4 fis. The 
detection efficiency for single photons depends directly on the fraction of 
the SPE charge distribution above trigger threshold. For high energy neu- 
trino events, many waveforms show contributions from multiple photons, all 
of which could provide useful information during event reconstruction. The 
overall light yield provides an estimate of the neutrino energy, and the space 
and time distribution of light helps to reconstruct direction and reject back- 
grounds. The time distribution of photons can be extracted from each PMT 
waveform if the response to single photons is well understood. The response 
to each photon is approximately given by the average SPE waveform, scaled 
randomly according to the complete charge probability distribution. 

In order to mimic the ambient temperature in the ice, PMTs were placed 
in a freezer box at —32 °C and illuminated by diffused light from a 375 nm 
UV LED. The light was generated in 10 ns pulses with intensity of about 
0.1 photons per shot (~0.02 photoelectrons per shot), dim enough to initiate 
only SPE signals. 

Figure E] shows the average SPE waveform, measured at the output of the 
AC coupling transformer with a digital storage oscilloscope (LeCroy LT374, 
500 MHz bandwidth, 0.5ns samples). Here the 95 Q input impedance of the 
DOM's front end amplifier was replaced by the series combination of a 50 Q 
resistor and the oscilloscope input. 
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Individual waveforms have different amplitudes but their shapes are simi- 
lar to within a few percent. The waveform is dominated by a peak of Gaussian 
shape (a = 3.2 ns) which accounts for 83% of the area. A tail on the late 
side of the peak accounts for the remaining area and exhibits a small amount 
of ringing. About 90% of the charge is collected before 10 ns after the peak. 
A substantial part of the observed pulse width is attributed to the damping 
resistors and the coupling transformer (Section [3]) . 
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Figure 3: Average of 10,000 SPE waveforms for one PMT at gain 1 x 10 7 , as seen at 
the secondary of the AC coupling transformer. Results for other PMT samples are very 
similar. The solid and dashed curves correspond to new and old transformer designs 
discussed in Section [SJ 

To study the total charge in SPE events, a computer-controlled integrat- 
ing ADC module (LeCroy 2249 A) was used to integrate charge in a 70 ns 
window, triggered by the synchronization signal of the LED pulse genera- 
tor. Figure H] shows a typical charge histogram, which exhibits a clear SPE 
peak to the right of the pedestal peak. The Gaussian part of the SPE peak 
corresponds to a charge resolution of approximately 30%. 

The non-Gaussian component rising below 0.3 times the SPE charge in 
Fig. H] has been studied to verify that such small pulses actually reflect in- 
time detection of photons, and not accidental coincidences of noise pulses 
such as from thermionic emission at the dynodes. The check for a noise 
contribution was done with the LED light output disabled (but not the syn- 
chronization signal that triggers acquisitions); all counts outside the narrow 
pedestal region were greatly suppressed compared to Fig. HI 

The low-charge component has been described in the past for many 



PMTs [2l|, and has been attributed to a sizable probability for backscat- 
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tering of the primary photoelectron at the first dynode 22|, |23| > leading to 
events where only a few secondaries are produced instead of the usual 10-20. 
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Figure 4: Typical pedestal-subtracted SPE charge histogram at gain 5 x 10 7 , including 
pedestal peak. Fits refer to Eq. (TTJ), with the constraint q T /qo = 0.2. The remaining 
parameters are optimized to fit this histogram for the curve labeled "Model" , while "Model 
(Scaled)" optimizes only the scale parameter qo while holding a q /qo and P e at values that 
describe the average of 120 PMTs. 



The shape of the low-charge component is important because even small 
pulses below the DOM's trigger threshold will be recorded in events with 
multiple photoelectrons. Therefore event reconstruction should account for 
the entire charge probability distribution down to zero charge, which we 



model Gaussian plus an exponential term [2jJ 
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Here P e is the fraction of events in the low-charge exponential part, go is 
the charge at the SPE peak which defines the PMT gain, <j q is the width of 
the Gaussian fit around the SPE peak, and q T is the decay constant in the 
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exponential component. Figure H] shows that this is a good model for the 
shape of the charge histogram away from the pedestal. 
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Figure 5: Scaled parameters from Eq. Q]as a function of PMT gain. The error bars show 
the la spread in parameters obtained for a sampling of 115 PMTs. 

Figure [5] shows results of fitting Eq. [I] for a large sample of PMTs at 
different gains above 5 x 10 7 , excluding the very low charge region q < 0.15go 
and the very high charge region more than 2a past the peak. The value 
of q T /qo is substantially degenerate with P e for describing observed spectra 
in the fitted range, so it has been fixed at the representative value of 0.20. 
The scaled quantities cr q /q , q T /qo, and P e are found not to vary strongly 
with the PMT gain. The very small pulses with q < 0.15go were omitted 
to avoid confusion with the tail of the pedestal distribution; results were 
the same if the low-charge cut was moved to 0.25go- The charge resolution 
cr q /q has been separately studied for gains between 10 7 and 10 8 and again 
no significant effects were seen. 

Figure [5] also shows the spread in parameters from PMT to PMT. The 
distribution in each parameter is approximately Gaussian, with the width 
shown by the error bars. The spread is substantial, but is not expected to 
have a large effect on data analysis, so the IceCube PMTs do not need to 
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be parameterized individually. Instead, an average model is currently used 
in event simulation and reconstruction, without modeling the spread. The 
similarity among PMTs can also be gauged from Fig. HI where data from one 
PMT is compared with a model curve scaled from the average fit results for 
120 PMTs. 

The above measurements were performed with diffuse light and represent 
an average over the photocathode surface. In a separate measurement at 
gain 10 7 , substantial differences were observed as a function of position; for 
example, the peak-to-valley ratio decreased to near unity close to the edge 
of the photocathode, compounding the effects of gain variation (Section [8]). 

7. Time resolution 

The timing of recorded SPE waveforms, relative to the photon arrival 
time, was studied at — 40 °C using fast pulses (FWHM 50 ps) from a Hama- 
matsu PLP-10 diode laser. Pulses were optically attenuated and diffused 
over the PMT face, yielding an average of 0.04 photoelectrons per shot. The 
wavelength was 405 nm. 

Each PMT was set for gain 10 7 based on its SPE charge spectrum. Hits 
greater than 0.4 times the SPE charge were recorded using the DOM digi- 
tization and readout electronics. Synchronization pulses from the laser were 
also digitized to indicate the true photon arrival times to within a fixed off- 
set. Hit times were defined as the points where each waveform reached 50% 
of its maximum, resulting in the time resolution histogram of Fig. [6J 

The main peak of the time histogram has width equivalent to a Gaussian 
of a = 2.0 ns, although the rising and falling edges of the peak fit better to 
half-Gaussians with a = 1.2 ns and 2.7 ns, respectively. Some of the width 
can be attributed to simultaneous illumination of the entire photocathode 
in our tests. When illuminated at the center only, the width decreased to 
1.5 ns; conversely, the outer 3 cm of the photocathode exhibited additional 
delay of about 3 ns and additional broadening. The data acquisition sys- 
tem contributed a time smearing of less than 0.6 ns, which has not been 
subtracted. 

About 4% of hits are found in a shoulder (25-65 ns) and secondary peak 
at 71ns, and 0.2% make up a corresponding tertiary structure (85-160 ns). 
The delayed hits are believed to arise when an electron trajectory is scattered 
back from the first dynode towards the photocathode, where it turns around 
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Figure 6: Typical distribution of SPE hit times for an IceCube PMT, when illuminated 
by narrow pulses from a diode laser. Counts are shown for 1 ns time bins, referenced to 
a laser synchronization pulse. A small fraction of late pulses (dashed line) are due to 
laser afterglow plus the random background count rate; hit fractions described in the text 
are corrected accordingly. The random background also explains the counts before the 
synchronization pulse. 



and then eventually arrives back at the first dynode to initiate the pulse 23 



Because of strong photon scattering in the ice, the dispersion of hit times 
by the PMT at the 2 ns scale is not a limiting factor for reconstruction in 
IceCube; likewise for the tail at late times. Considering the spacing between 
DOMs, photons must typically travel tens of meters before detection, which is 
comparable to the scattering length of around 25 m Q. A detailed simulation 
of photon scattering 2(| showed that at 10 m distance, about 40% of photons 
are delayed by more than 5 ns, and 10% of photons were delayed between 
20 ns and 80 ns. This is larger than the corresponding effects from the PMT 
itself. On the other hand, DOMs close to a high energy track can be expected 
to detect at least one photon with negligible delay, and then the very small 
1.2 ns dispersion on the early side of the time resolution peak may be relevant 
when reconstructing arrival time of the earliest photon or the pulse rise time. 

The time resolution study also reveals DOM-to-DOM differences in the 
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nominal delay of SPE waveforms relative to photon arrival time. This delay 
includes PMT transit time plus signal delays between the PMT output and 
the digitizer. The PMT transit time is found to vary according to the square 
root of the applied voltage, 



where k = 0.017 ns/volt is the slope at V = 1500 V. The voltage applied 
to each PMT is set for a design gain 10 7 , which is achieved between 1050 V 
and 1600 V in 99.9% of PMTs. The resulting RMS spread of the overall time 
offset is 2.7ns. We find 5% of DOMs more than 5ns away from the mean, 
so the DOM-to-DOM corrections are currently included in reconstruction. 

8. Two-dimensional photocathode scan 

The number of photons arriving at the PMT is calculated from the ob- 
served photoelectron signals via the PMT optical efficiency. This is separated 
into an overall "absolute efficiency" and an angular dependence. The domi- 
nant factor in angular dependence is just the amount of photocathode area 
which can be seen from various directions. However this has to be adjusted 
for the fact that the photocathode surface is very large and different portions 
do not all yield the same efficiency. 

We have systematically analyzed the variation of efficiency with photo- 
cathode position at 25 °C, using a two-dimensional scan system. A UV LED 
(370 nm) with collimator produced a 1 mm spot which was moved along the 
curved PMT surface, maintaining normal incidence of the light. The LED 
delivered approximately 125 photons per 80 ns pulse. The anode voltage was 
set for gain 10 7 at the center of the photocathode, as measured by the SPE 
charge peak (q in Eq. [TJ). The PMT pulse charge for each position was then 
measured by an integrating ADC triggered by the LED pulser. 

Figure [7] shows typical response maps on the cathode surface. The mea- 
sured charge is proportional to the net photomultiplier efficiency and reflects 
the combined position dependence of photocathode quantum efficiency, col- 
lection efficiency, and dynode multiplication. PMT to PMT variation of the 
efficiency at a given spot on the photocathode may be as great as 40%, how- 
ever, the spread in the area-integrated efficiency from PMT to PMT is much 
smaller, of order 10% (see Sec. |9]). The average map shows a uniform falloff 
in the edge region, except for a small bias in the +X direction. All PMTs 
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Figure 7: Position dependence of the light pulse response for three example PMTs, and 
the average of 135 PMTs (lower right). The X-Y coordinates measure distance from the 
center of the photocathode along the curved PMT face. The color scale indicates the PMT 
output pulse charge in units of the SPE charge, averaged over many pulses. 
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were measured in the same orientation, so this bias could be associated to 
the first dynode position or the geomagnetic field. 

A small part of the variation seen in the scans can be attributed to sys- 
tematic errors, which arise mainly from the geomagnetic field and LED lumi- 
nosity variance. The geomagnetic field of 462 mG is attenuated by about 50% 
with a shield made of /x-metal sheet and wire, as also used in IceCube DOMs. 
By comparing measurements with the PMT rotated from its standard ori- 
entation in various ways, we determined the field's effect on the efficiency 
variation is about 10%. Because the magnetic shield is the same, the overall 
results are representative of what is expected for deployed IceCube DOMs. 
(The field at the South Pole is 553 mG, and more vertical relative to the 
PMT axis.) The time dependence of LED luminosity affected the shape of 
each scan by less than 2%, as seen by reproducibility of the scan results. 

By reducing the light intensity to give only SPE hits, a similar map has 
been constructed for gain variation. The gain can vary as a function of po- 
sition because the corresponding photoelectron trajectories arrive differently 
at the first dynode, leading to different yields of secondary electrons. The 
observed gain (defined by go m Eq. [T]) varies within ±10% over the active 
region when high voltage is set for gain 5 x 10 7 at the center. However, 
the low-charge contribution to the SPE charge spectrum (P e in Eq. [1]) was 
found to also vary with position, so that the peak-to-valley ratio decreases 
close to unity near the edge of the photocathode. In this way, the average 
charge delivered per photoelectron was observed to decrease by up to 30% at 
nominal gain 10 7 . Because of these effects, the detection efficiency map for 
single photons using a specific discriminator threshold can differ somewhat 
from the maps of Fig. [7J 

The integrated sensitivity for a broad beam of photons incident from 
a particular direction follows from the cathode efficiency maps by averag- 
ing over the surface seen from that direction [27]] , with correction for non- 
normal incidence on surface elements as appropriate [28j]. For this purpose 
the relative efficiency map is assumed not to vary with wavelength, i.e., each 
position is assumed to obey the spectral response curve given by the man- 



ufacturer [16j. The averaging substantially reduces the effect of variations 
over the surface so the sensitivity is not strongly dependent on direction at 
moderate polar angles. Only light that arrives at large polar angles relative 
to the PMT axis will primarily illuminate the equator region and therefore 
show strong azimuthal dependence. Likewise, the variation in charge spec- 
trum from center to edge has little effect after averaging. 
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The scans were performed on a small fraction of the IceCube PMTs, so 
only the average variation with polar angle is used in simulation and recon- 
struction. For IceCube, the remaining PMT-to-PMT variation in directional 
sensitivity has very small consequences, because light is typically scattered 
after traveling about 25 m through the ice, and additional scattering takes 
place in the refrozen ice in the hole where the DOMs are deployed. The 
PMT-to-PMT variation, as well as the position dependence itself, could be 
more important for detectors deployed in water where scattering lengths are 



much longer [29 



9. Absolute efficiency calibration 

The absolute calibration of PMT optical efficiency is important because 
IceCube uses the observed number of photons to estimate energy in recon- 
structed neutrino interactions. Showers initiated by electrons or tau leptons 
yield light in proportion to the energy, and so do muons above 1 TeV where 
energy loss is dominated by direct pair production, photonuclear interactions 
and bremsstrahlung 0, 0] . 

Optical efficiency can be studied after deployment by using light from 
muons (produced in cosmic ray showers above IceCube) or from calibrated 
beacons deployed in the ice nearby. However it is hard to isolate the PMT 
response from the effects of light scattering and attenuation in the ice, which 
have some uncertainties [?|]. 

Here we describe the laboratory calibration of standard PMTs installed in 
16 IceCube DOMs distributed throughout the array. The calibrated PMTs 
provide direct information for the energy calibration of IceCube, and will also 
help clarify the ice effects in other studies, which is an important subject on 
its own. 

9.1. Technique 

Our setup for measuring a PMT's UV photon detection efficiency is shown 
in Fig.[HJ A pulsed 337 nm laser beam is passed through a chamber containing 
pure nitrogen gas, and the PMT to be calibrated is illuminated by the tiny 
amount of light that is Rayleigh scattered at about 90 degrees from the 
beam. The PMT is rotated inside the dark box to probe different positions 
on the photocathode surface. The primary beam intensity is measured with 
a calibrated silicon photodiode "energy probe", and sets the fundamental 
scale for our efficiency measurement. A pressure gauge and temperature 
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sensor provide the gas density. Beam intensity, geometry and gas density are 
then folded with the well-known Rayleigh scattering cross section to obtain 
the absolute number of photons per pulse incident on the PMT. Individual 
photon detections are counted in each pulse and divided by the number 
incident to obtain the optical sensitivity at 337 nm. Additional corrections 
(Sec. I9.4p are needed to obtain DOM efficiency at wavelengths around 400 nm 
where IceCube is most sensitive. 

The measurement combines effects of quantum efficiency and collection 
efficiency and can be directly applied in IceCube analysis. It is different 
from the usual quantum efficiency measurement, which is based on cathode 
current response to a calibrated DC light source 



The laser (Spectra-Physics VSL-337ND-S) emits 4 ns pulses containing 
~10 10 photons, as measured by the silicon photodiode probe (Laser Probe, 
Inc., RjP-465). After a warm-up delay, pulse energies are stable to within 
±2%. ' 

The beam width is 1 mm, so the illuminated gas volume may be consid- 
ered as a line source of Rayleigh scattered light. Apertures between the beam 
and the PMT define a source region with effective length of about 1 cm and 
a spot size on the PMT of about 1.5 cm. Photons can reach the PMT if they 
are scattered from the source region into a solid angle of about 7.6 x 10~ 4 sr 
around the 90° direction, while other scattered photons are eventually ab- 
sorbed on baffles or other surfaces inside the chamber. 

The Rayleigh scattering cross section for a circularly polarized beam on 



nitrogen gas is taken as [31 



" (1 + cos 2 6) x (3.50 ± 0.02) x lO" 26 cm 2 (3) 



dVt 16tt 



with 9 as the polar angle relative to the beam direction. The geometrical 
integration over the source region and corresponding solid angles is handled 
in a detailed ray-tracing calculation. After accounting for pressure and tem- 
perature, this yields the overall number of scattered photons reaching the 
PMT, typically 0.5 per pulse. With detection efficiency around 20%, this 
corresponds to ~0.1 SPE per pulse. 

For counting photon detections, the PMT output charge is integrated for 
each laser pulse with a CAMAC ADC. The gating time is 184 ns which is long 
enough to include the late PMT pulses described in Section [71 The PMT gain 
is set close to 10 8 as defined by the SPE peak, q in Eq. [TJ We then count the 
number of events with charge q greater than a threshold q t h = 0.5go, which 
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Figure 8: Schematic view of the absolute calibration system. Photodiode 1 establishes the 
beam intensity, which is used to predict the amount of Rayleigh-scattered light reaching 
the R7081-02 PMT. The Hamamatsu H7195 PMTs and the other photodiode are for 
monitoring. Optical baffles define the scattering geometry, and chamber windows are 
anti-reflection coated UV quartz (R < 0.1% at 337 nm). 
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Table 2: Measured photon detection efficiency (77) and photon effective area (A c g) at 
25 °C for four different PMTs at wavelength 337 nm and gain 10 8 . Values for g t h = were 
extrapolated using Eq. [T] where model parameters were fit independently for each PMT. 
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can be clearly discriminated in the charge histogram. A small correction for 
events with multiple photoelectrons yields the number of detected photons 
with q > q th . 

The PMT efficiency t] for g t h = 0.5go is then given by comparing the 
number of detected photons to the number reaching the PMT. The efficiency 
for other charge thresholds can be computed by extrapolation with the SPE 
charge response model, Eq. [TJ 

9.2. Results 

Figure IH] shows the measured detection efficiency as a function of distance 
from the cathode center. It also shows that the absolute efficiency measure- 
ments follow closely the shape expected from the 2D relative efficiency scans 
(Section [8]). Consequently, the relative efficiency scans can be normalized 
to the absolute measurements and used to estimate the absolute efficiency 
averaged over any given area of the photocathode surface. 

Table [2] lists measured efficiencies at the center and averaged over the 
whole photocathode area. The latter is defined to include all points within 
15 cm of the PMT axis, measured along the curved surface. Table [2] also 
includes the photon "effective area" which means the amount of ideal surface 
with 100% efficiency that corresponds to the actual convolution of area and 
efficiency. Here it is quoted for light uniformly spread over the PMT surface, 
with normal incidence. For use in IceCube analysis, a similar calculation is 
performed for unidirectional beams as a function of the beam angle, folding in 
the variation of cathode response and optical effects at material boundaries. 

The detection efficiencies in the central area of the photocathode are 
close to 20% if extrapolated to q t h = 0. These values have been compared 
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Figure 9: Detection efficiency as a function of distance from the cathode center for three 
different PMTs. The points with error bars show absolute efficiency measurements ex- 
trapolated to qth — 0. The histogram curves show corresponding results from the 2D 
relative efficiency scans (Section [8]) after normalization. Systematic uncertainties for the 
histograms are described in Section [5J 



on a PMT-by-PMT basis with measurements by Hamamatsu using cathode 
response to DC light sources. We find very good agreement, which implies 
that the collection efficiency is not much less than 100% at the cathode center. 

IceCube PMTs operate at lower voltages than used in this measurement 
(gain 10 7 instead of 10 8 ), so collection efficiency is expected to be slightly 
lower. However, this effect is expected to be concentrated at the edges where 
the electron optics are less ideal, and the efficiency falloff near the edges is 
obtained from the 2D relative efficiency scans (Section [8]). Since those scans 
were done at gain 10 7 , no additional correction is necessary for the results in 
Table H 

9.3. Uncertainties 

The overall systematic uncertainty Arj/r] of the PMT detection efficiency 
measurement is 7.7%, as detailed in Table |3j In addition, the measurement 
of each position on the PMT face has a typical statistical uncertainty of 
about 5%, set by the number of SPE hits recorded. This is reduced to about 
2% when calculating efficiency for the whole PMT by combining information 
from the individual face positions, but the extrapolation relies on the 2D 
map (Section [8]) which has comparable uncertainties. 

The dominant contributions to systematic error arise from the laser beam 
energy measurement, the geometry of the scattering chamber, and the geo- 
magnetic field. The first two enter directly into the calculation of the number 
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Table 3: Systematic error budget for the PMT efficiency calibration. 



Source 


Arj/i] 


Laser beam energy 


5 % 


Aperture 


4 % 


Ambient magnetic field 


4% 


Pressure and temperature 


1 % 


Polarization 


1 % 


Rayleigh cross section 


0.5 % 


Dark noise / cosmic rays 


0.2 % 


Overall 


7.7% 



of photons reaching the PMT. The beam energy comes from the laser energy 
probe, which was factory calibrated to within 5% at 337 nm. The aperture 
uncertainty of 4% comes from geometrical survey of the chamber, which is 
used in the ray tracing program. 

The ambient geomagnetic field is 462 mG and is unshielded in the current 
setup. By changing the orientation of the PMT, we showed that it affects 
the point-to-point response map at the 20% level but the average over the 
surface varies by only 4%. 

The Rayleigh scattering angular distribution depends on the polarization 



of the laser beam [32J, so care was required to limit any polarization effect. 
The effect can be strong because the PMT only sees a Rayleigh scattering 
signal from the vertical component of polarization: the horizontal compo- 
nent induces a dipole moment oscillating perpendicular to the beam in the 
horizontal plane, which cannot emit power to the PMT which is in the same 
direction. On the other hand the energy probe reads the total power regard- 
less of polarization, so the fraction of power in the vertical direction must 
be under control. In our setup (Fig. [S]), the laser beam is first linearly po- 
larized at 45° and then passed through a quartz A/4 waveplate to convert 
it to 100% circular polarization. By rotating a linear polarizer in the beam, 
we verified that the resulting horizontal and vertical components are equal 
to within better than 1%, which leads to a limit of 1% for the corresponding 
systematic effect on the efficiency measurement. 

Several analyses were performed to show that the scattered photons are 
coming from Rayleigh scattering and not other sources such as residual sus- 



27 



pended dust. We checked for expected scaling with gas density down to low 
pressure, symmetry of scattered light seen in forward and backward moni- 
toring PMTs (see Fig. [8]), and repeatability of measurements after long time 
intervals. 

The main systematic uncertainties in our method could be reduced if 
desired, comparable to the current statistical precision of about 2%. To 
accomplish this, one would calibrate the silicon photodiode "energy meter" 
at the 1% level; measure the aperture geometry more precisely; and provide 
good shielding from the geomagnetic field. 

9-4- Additional corrections 

There are additional steps to obtain the detection efficiency of DOMs 
from the PMT efficiency measurements, and these will be reported sepa- 
rately along with other studies on assembled DOMs. Corrections include 
attenuation of light in the glass pressure housing and the gel used for optical 
and mechanical coupling, wavelen gth dependence in both the PMT sensitiv- 



ity (quoted by the manufacturer |l6|) and the attenuation factors, and the 



geometry of incident rays. These effects are included in a detailed optical 



simulation of the DOMs [27], |33| which will be compared to laboratory mea- 
surements on assembled DOMs. A full detector simulation can be used to 
combine the absolute efficiencies at 337 nm, the wavelength dependences, 
and the spectrum of light received from neutrino interactions. Investigations 
of the combined effect show that IceCube detects signal photons in a broad 
range centered on about 400 nm. These studies will be presented elsewhere. 

Our measurements were at 25 °C, and some temperature dependence 
can be expected. The manufacturer quotes a temperature coefficient of 
— 0.2%/°C for cathode sensitivity (l6| . This is being directly addressed by 
relative measurements of optical efficiency of assembled DOMs at — 45 °C, 
—20 °C and room temperature. 



10. PMT linearity and saturation behavior 

For most neutrino interactions expected in IceCube, any one PMT will 
not detect more than a handful of photons. For such events, and even when 
up to a few hundred are detected, the PMT output is proportional to the 
number of photons detected. However, some of the most interesting signal 
events would be expected to deposit large amounts of energy within tens of 
meters of individual PMTs, and then the PMT response can be less than 
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10 1 10 2 10 3 10 4 

I = Ideal Current (mA) 

Figure 10: (a) PMT saturation curve for gain 10 7 . The measured current is plotted against 
the instantaneous light level, defined in terms of the current that would be expected for 
an ideal (linear) device. Data points correspond to peak currents for the indicated pulse 
widths. The fit curve is given by Eq. |4j with parameters chosen optimally for this PMT 
(serial number SA2747). (b) Effect of PMT gain on the saturation curve. Data points 
correspond to peak currents for 200 ns pulses for PMT serial number AA0020. The curves 
show the fitted parameterization, Eq. [5] 
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proportional. Optimal reconstruction requires measuring the linearity limit 
and modeling the nonlinear saturation behavior. 

To study saturation behavior, PMTs were illuminated with LED pulses 
of various durations and brightnesses. The duration of the light pulses was 
varied from 3ns to 1 /is, as measured with a fast PIN photodiode (Ins re- 
sponse time). A set of calibrated neutral density filters was used to control 
the light level. For a given LED brightness, illumination was first measured 
by the PMT signal, using a filter with sufficient attenuation to allow linear 
operation of the PMT. The observed signal was converted to a photoelectron 
rate and a total number of photoelectrons (p.e.) using the SPE charge go, 
determined in a separate step. Then the illumination level was increased by 
using different filters, with the new number of p.e. calculated from the ratio 
of filter attenuation coefficients. 

Figure ITOf a) shows the observed peak anode current / as a function of 
the ideal peak current Io, defined as the peak p.e. rate times the SPE charge. 
At gain 10 7 , the PMT response is linear within 10% up to currents of about 
50mA (31p.e./ns), but saturates completely at about 150mA. Peak responses 
to different light pulse widths from 3 ns to 400 ns lie along a single curve. The 
3 ns and 20 ns width pulses were approximately Gaussian in shape, so the 
observation of identical peak response supports a saturation model where the 
observed current is a direct function of the instantaneous illumination, with 
little cumulative effect from previous illumination. In particular the data 
are inconsistent with models expressed in terms of total pulse charge, which 
were used in some older versions of the IceCube simulation software. The 
400 ns pulses were approximately rectangular in shape and the output current 
mirrored this shape well even in the saturation region, again as expected 
for an instantaneous current saturation model. Even long light pulses near 
saturation level show only about 5% decline from 100 mA after 1 fis of steady 
illumination. (Note this small history effect is independent of the transient 
gain loss due to discharge of the dynode capacitors, which remains below 1% 
for such a pulse.) 

The same saturation behavior was found to apply regardless of what 
part of the cathode was illuminated, even at —30 °C, which indicates that 



photocathode surface resistance [16( is not important on the relevant time 
scales. 

The instantaneous current response is well parameterized by the follow- 
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Table 4: Saturation curve parameters for three PMT samples, as defined for Eq. [4] 



PMT Serial No. 


A (rriA) 


B 


C 


AA0020 


126 


2.02 


2.98 


SA2747 


138 


2.05 


3.23 


SA2749 


138 


1.82 


2.67 



rng: 

The parameters A, B and C differ substantially from one PMT to another 
(Table Hj), so the model should not be used to invert observed saturated 
pulses unless each PMT is fully characterized. 

Figure [TOT b) shows additional measurements at a range of lower gains 
down to 10 5 , relevant for IceTop DOMs. The model of Eq. H] continues to 
apply over the full range if the parameters are scaled approximately as powers 
of the gain, as shown by the curves which are scaled by 7 = Gain/10 6 : 

A{i) = 285 7 a5 7(l + 7 1/4 ) 2 
B(j) = 13 7 °- 18 /(1 + 7 1/4 ) 2 

C( 7 ) = 0.32 7 - ai3 (l + 7 1 / 4 ) 2 (5) 

The given parameters apply to a single measured PMT, but similar scaling 
behavior can be expected for other examples; as a first estimate one would 
adjust the leading coefficients in each parameter to match measurements at 
a particular gain, and retain the same scaling with gain. Note the good 
numerical behavior of the scaling equations allows them to be used also for 
estimates outside the given range of gain. 

The instantaneous current model also helps understand how the response 
to narrow light pulses (3 ns FWHM) broadens as intensity increases (Fig. fTTT) . 
With the light pulse strongly attenuated (220 p. e.), the PMT output pulse 
width is similar to the SPE response, about 10 ns. As more light is allowed 
to reach the PMT (3,700 p. e.), first a gradual broadening occurs to about 
20 ns width. This broadening follows from Eq. H] because the peak current 
is more attenuated than the rising and trailing edges. At this point a tail 
is visible, along with a second peak delayed by about 60 ns relative to the 
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Figure 11: Average waveforms observed in PMT serial number SA2747 for 3 ns (FWHM) 
light pulses with progressively higher intensity: (a) main peak; (b) secondary peak due to 
unusual electron trajectories; (c) pre-pulse. 
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main peak. These are consistent with the late photoelectron responses seen 
in SPE time resolution measurements (Section [7]), except that the relative 
sizes of main peak and tail are altered by saturation in the main peak. At 
still higher light levels (210,000 p. e.), the second tail peak is comparable in 
size to the fully saturated main peak, and the total width is dominated by 
the combination of the two peaks. 

The highest light level in Fig.[TT]also exposes a small pre-pulse 30 ns before 
the main peak, as well as a substantial afterpulse starting several hundred 
ns later (see Section EJ. The pre-pulse is ascribed to photoelectrons ejected 
from the first dynode, and is somewhat exaggerated in Fig. [11] because the 
light source was aimed at the center of the cathode with the dynode directly 
behind. The individual quanta comprising the pre-pulse were separately 
studied using SPE-level illumination, and were found to be between 1/10 and 
1/20 of the SPE pulse size, occurring at less than 1% of the SPE rate. The 
ratio between sizes of SPE pulses and pre-pulse quanta is similar to the typical 
first dynode gain reported by the manufacturer 24j. Because individual pre- 
pulse quanta are below threshold for triggering DOMs, they have a small 
impact on event reconstruction. The combined pre-pulses from many photons 
would only be observable for a large, narrow light pulse ( > 5000 photons 
detected within 30 ns). Pulses originating more than ~25m from a DOM 
would generally be broader than this, due to scattering in the ice 0. 

The saturation model (Eq. H]) can be important for reconstruction of 
very high energy neutrinos that produce electromagnetic or hadronic show- 
ers. Ideally, reconstruction would rely most heavily on the PMTs closest to 
a shower, because the light pulse is broadened and attenuated as it travels 
through the ice Q; however these PMTs can be saturated for high energy 
events. The energy where saturation effects become important can be esti- 
mated by choosing a characteristic distance of 60m, which is about half the 



inter-string spacing. At this distance, simulations [20] show that a 600 TeV 
shower yields peak intensity of 30p.e./ns, equivalent to the linearity limit of 
50 niA. Above this energy, signals in close PMTs require a correction for sat- 
uration. Above ~10PeV, many nearby PMTs are badly saturated and the 
shower energy measurement must rely mostly on far-away PMTs. However, 
even badly saturated PMTs measure the beginning and end of the pulse, 
which can be used to constrain the event geometry. 



33 



0.0 



-2.0 - 



-4.0 - 



g - 6 - 



Z3 

o 

CD 



Z3 
CO 

CD ru 



-80 - 



-1 20 



-, 




3.0 x 10 3 PEs 
1.6 x 10 4 PEs 

1.1 x 10 5 PEs 
—i 1 1 — 





2.0 x 10 5 PEs - 
1.0 x 10 6 PEs 
4.4 x 10 6 PEs 



10 



Figure 12: PMT afterpulse waveforms for bright flashes lasting 40 ns. Each curve is 
averaged over many flashes. The primary response (a) is off-scale and saturated at most 
of these intensities. Brightness of flashes is determined by the method of Section[Tni which 
is independent of PMT saturation effects. Prominent afterpulse peaks are seen around (b) 
600 ns, (c) 2 fis and (d) 8 /js after the primary response. The afterpulse waveforms grow 
nearly linearly with flash intensity up to the maximum intensity measured. The PMT was 
operated at gain of 10 7 , with the anode at 1326 V. Some variation in afterpulse waveforms 
is observed from one PMT to another of the same type. Measurement of tails at very late 
times is slightly affected by an AC coupling time constant of 80 fxs. 
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11. Afterpulses 



As shown in Figs. [6] and [TT], the prompt response to a light pulse has a 
tail extending to about 100 ns. Afterpulses are seen in the range of 300 ns to 
11 /is. Such afterpulses are a common feature of PMTs, and are attributed 
to ionization of residual gases by electrons accelerated in the space between 



dynodes [34J . Ions created in this way can be accelerated back to the photo- 
cathode, causing ejection of electrons which are subsequently amplified like 
the original photoelectrons. Some ions strike one of the dynodes instead, 
but the corresponding ejecta are amplified much less and could easily go 
undetected. 

Afterpulse measurements were made at 25 °C with LED pulses of 40 ns 
width, using calibrated optical attenuators to control the intensity, as for the 
saturation measurements (Section [TO)) . In a bright LED flash, many individ- 
ual ions are created, and their responses add up to an afterpulse waveform 
with well defined peaks and valleys (Fig. IT2"j) . The various peaks are believed 
to correspond to ions of different masses, according to their individual flight 
times in the accelerating field 34j. Prominent afterpulse peaks for this PMT 
occur around 600 ns, 2 /xs and 8 /xs after the main response peak. The peaks 
are fairly wide and no period is entirely devoid of afterpulses until after 11 /xs. 

The average afterpulse waveform grows almost linearly with the flash 
brightness even up to the highest intensity studied (4.4 x 10 6 p.e. in 40 ns), 
where the primary response is completely saturated at ~1000p.e./ns (Fig. [TO]) . 
This suggests that observed afterpulses arise primarily from ions generated 
in earlier stages of the multiplier, whose electron currents continue to rise 
even when later stages have saturated. 

Up to primary pulses of 1 x 10 6 p.e., the integral from 300 ns to 11 /xs 
corresponds to 0.06 SPE per primary photoelectron. 

For dimmer flashes, individual events have a small number of afterpulse 
electrons. These appear as separate single afterpulses distributed in time, 
with probability that can be approximately predicted from the average wave- 
forms of Fig. d2J Because different ions are associated with different time 
ranges, and because some ions eject multiple electrons from the photocath- 
ode, each afterpulse delay range will be characterized by a different funda- 
mental charge distribution. We have observed corresponding peak charges 
from 1 SPE to 13 SPE, consistent with a recent more detailed study of indi- 
vidual ion afterpulses [iil ]. 

The above observations are from study of only a few PMT samples, and 
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the numbers quoted pertain to only one (serial AA0020). Although quantita- 
tive differences are seen from one sample to another, the information allows 
one to assess whether afterpulses affect IceCube event reconstructions, and 
to limit small systematic errors. If a particular physics analysis then appears 
sensitive to afterpulses, a larger sample of PMTs would have to be studied 
quantitatively to provide the necessary corrections. 

Given the small ratio of charge between afterpulse and primary pulse, it 
can be expected that most IceCube analyses will not be strongly sensitive 
to the details of afterpulses. Typical IceCube events yield hits in each PMT 
that are spread over times of a few hundred nanoseconds, well before the 
main part of the afterpulse distribution. For very high energy events (e.g. 
electron energy 1 PeV in the deep ice), signals are likely to be seen by PMTs 
500 m away where arrival times are dispersed over 2/xs (FWHM), and then 
the afterpulse distribution becomes more relevant. However, the main effect 
is a minor distortion of the late part of the pulse, which already has an 
intrinsically long tail due to scattering. Generally one does not lose much 
information by disregarding details of the waveforms at late times. 

However, some events can have multiple peaks in the photon time dis- 
tributions, and then a characterization of afterpulses can be important for 
proper reconstruction. The most common case is an event with coincident 
arrival of one or more downgoing muons from cosmic ray showers above 
the detector, which calls for disentangling the hits originating from multiple 
tracks, and therefore also the afterpulses. Multiple muons can also arise from 
a single shower, and when the resulting tracks are well separated they can 
yield multiple hits. More intriguing is the possibility of v T interactions which 
can create two showers of particles separated by hundreds of meters (36| . In 
such a case some PMTs can see pulses of light separated by a few microsec- 
onds, so effects of afterpulses should be considered carefully. The late pulses 
described in Section [7] should also be considered in these contexts. 

12. Summary 

The R7081-02 PMT has been characterized and key findings were dis- 
cussed in the context of IceCube physics goals. We observe a single-photoelectron 
time resolution of 2.0 ns averaged over the face of the PMT. A small fraction 
of the pulses arrive much later, with about 4% between 25 and 65 nsec late. 
We also observe prepulsing and afterpulsing, with afterpulsing occurring up 
to 11 [is late. The single photoelectron charge spectrum is well fit by a Gaus- 
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sian corresponding to charge resolution near 30%, plus a contribution at low 
charge which is represented by an exponential. The dark rate was measured 
to be 300 Hz in the temperature range —40 °C to —20 °C. A new method for 
optical sensitivity calibration has been demonstrated, which uses Rayleigh 
scattering to scale from the intensity of a primary laser beam to the much 
smaller number of photons reaching a target PMT. Measurements of dark 
rate, single photon detection efficiency, single photoelectron waveform and 
charge, time resolution, large pulse response, and afterpulses will serve as 
input for detailed simulation of IceCube physics events. 



Acknowledgment s 

We acknowledge support from the following agencies: U.S. National Sci- 
ence Foundation-Office of Polar Program, U.S. National Science Foundation- 
Physics Division, University of Wisconsin Alumni Research Foundation, U.S. 
Department of Energy, and National Energy Research Scientific Computing 
Center, the Louisiana Optical Network Initiative (LONI) grid computing 
resources; Swedish Research Council, Swedish Polar Research Secretariat, 
and Knut and Alice Wallenberg Foundation, Sweden; German Ministry for 
Education and Research (BMBF), Deutsche Forschungsgemeinschaft (DFG), 
Research Department of Plasmas with Complex Interactions (Bochum), Ger- 
many; Fund for Scientific Research (FNRS-FWO), FWO Odysseus programme, 
Flanders Institute to encourage scientific and technological research in indus- 
try (IWT), Belgian Federal Science Policy Office (Belspo); Marsden Fund, 
New Zealand; Japan Society for the Promotion of Science (JSPS); M. Ribordy 
acknowledges the support of the SNF (Switzerland); A. Kappes and A. GroB 
acknowledge support by the EU Marie Curie OIF Program; J. P. Rodrigues 
acknowledge support by the Capes Foundation, Ministry of Education of 
Brazil. 



References 

[1] IceCube Collaboration: J. Ahrens et al, Astropart. Phys. 20 (2004) 507 
[2] IceCube Collaboration: J. Ahrens et al, "IceCube Preliminary Design 



Document" (2001) http://icecube.wisc.edu/ 



[3] T.K. Gaisser, F. Halzen, and T. Stanev, Phys. Rep. 258 (1995) 173; 
J.G. Learned and K. Mannheim, Annu. Rev. Nucl. Part. Sci. 50 (2000) 
679; F. Halzen and D. Hooper, Rep. Prog. Phys. 65 (2002) 1025 



37 



[4] A. Karle for the IceCube Collaboration, "IceCube", Proceedings of the 
31st International Cosmic Ray Conference, 7-15 July 2009, Lodz, Poland 

[5] W. Lohmann, R. Kopp and R. Voss, "Energy Loss of Muons in the 
Energy Range 1-10000 GeV", CERN Report 85-03 (1985) 



[6] D. Chirkin and W. Rhode, preprint arXiv: hep-ph/040 70 75 v2 (2008) 



[7] Ackermann, M., et ai, J. Geophys. Res., 11 (2006) D13203 

D. Chirkin for the IceCube Collaboration, Proceedings of the 30th In- 
ternational Cosmic Ray Conference, 3-11 July 2007, Merida, Yucatan, 
Mexico, session HE1.5; in arXiv:0711.0353 [astro-ph], 151-154 

[9] AMANDA Collaboration: E. Andres et a/., Astropart. Phys. 13 (2000) 
1 

[10] AMANDA Collaboration: J. Ahrens et a/., Phys. Rev. D 66 (2002) 12005 

[11] E. Middell, J. McCartin and M. D'Agostino for the IceCube Collabora- 
tion, Proceedings of the 31st International Cosmic Ray Conference, 7-15 
July 2009, Lodz, Poland 

[12] IceCube Collaboration: T. Gaiser et ai, "Ice Top Preliminary Design 



Document" (2003) http://icecube.wisc.edu/ 



[13] IceCube Collaboration: A. Achterberg et al, Astropart. Phys. 26 (2006) 
155 

[14] IceCube Collaboration: R. Abbasi et ai, Nucl. Instrum. Methods Phys. 
Res. A 601 (2009) 294 

[15] Hamamatsu Corp., R7081-02 Photomultiplier Tube Data (2003) 

[16] Hamamatsu Corp., Photomultiplier Tubes: Basics and Applications, 3rd 
ed. (2007) 

[17] Ankowski, A. et ai, Nucl. Instrum. Methods Phys. Res. A 556 (2006) 
146; Nikkei, J.A., Lippincott, W.H., and McKinsey, D.N., J. Instrum. 2 
(2007) 1004; Meyer, H.O., preprint arXiv:0805:0771vl [nucl-ex] (2008) 

[18] A. S. Dighe, M. Th. Keil and G. G. Raffelt, J. Cosmol. and Astropart. 
Phys. 0306 (2003) 005 



38 



[19] EMCO High Voltage Corporation, Sutter Creek, CA 

[20] Lundberg, J., et al, Nucl. Instrum. Methods Phys. Res. A 581 (2007) 
619 

[21] Dossi, R. et al, Nucl. Instrum. Methods Phys. Res. A 451 (2000) 623 

[22] Coates, P. B., J. Phys. D: Appl. Phys. 6 (1973) 153 

[23] Smirnov, O. Ju., Lombardi, P., and Ranucci, C, Instrum. Exp. Tech. 
47 (2004) 69 

[24] Hamamatsu Corp., Personal communication 

[25] Lubsandorzhiev, B.K., Pokhil, P.G., Vasiljev, R.V., and Wright, A.G., 
Nucl. Instrum. Methods Phys. Res. A 442 (2000) 452 

[26] A. Karle, "Monte Carlo Simulation of Photon Transport and Detection 
in Deep Ice: Muons and Cascades," in Proceedings of the Workshop on 
Simulation and Analysis Methods for Large Neutrino Telescopes, DESY 
1998 



[27] K. Hoshina, http://www.ppl.phys.chiba-u.jp/ROMEO/ 



[28] Motta, D., and Schonert, S., Nucl. Instrum. Methods Phys. Res. A 539 
(2005) 217 

[29] Fukuda, S., et al (The Super-Kamiokande Collaboration), Nucl. In- 
strum. Methods Phys. Res. A 501 (2003) 418; Aguilar, A., et al, As- 
tropart. Phys. 23 (2005) 131 

[30] R. Mirzoyan et al, Nucl. Instrum. Methods Phys. Res. A 572 (2007) 449 

[31] H. Naus and W. Ubachs, Opt. Lett. 25 (2000) 347-349, and references 
therein 

[32] W. Reckers, Y. Gu, E. W. Rothe and H. Voges, Appl. Spectrosc. 51 
(1997) 1012 

[33] S. Agostinelli, et al, Nucl. Instrum. Methods Phys. Res. A 506 (2003) 
250-303. 

[34] P. B. Coates, J. Phys. D 6 (1973) 1159-1166. 



39 



[35] K. J. Ma, et al, preprint |aiXivi091L5336vl [physics.ins-det] (2009) 
[36] J. G. Learned and S. Pakvasa, Astropart. Phys. 3 (1995) 267 



40 



